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We study the propagation of cold-atom wave packets in an interferometer with a Mach-Zehnder 
topology based on the dynamical phase of Bloch oscillation in a weakly forced optical lattice with a 
narrow potential barrier that functions as a cold-atom wave packet splitter. We calculate analytically 
the atomic wave function, and show that the expected number of atoms in the two outputs of the 
interferometer oscillates rapidly as a function of the angle between the potential barrier and the 
forcing direction with period proportional to the external potential difference across a lattice spacing 
divided by the lattice band energy scale. The interferometer can be used as a high precision force 
probe whose principle of operation is different from current interferometers based on the overall 
position of Bloch oscillating wave packets. 



I. INTRODUCTION 

The dynamics of a quantum particle in a periodic envi- 
ronment is in many ways similar to that of a free particle, 
but when subject to a uniform force it does not accelerate 
indefinitely but performs periodic motion — Bloch oscilla- 
tions. This effect is often explained using interference of 
Wannicr-Stark resonances, but can also be understood 
in terms of classical dynamics generated by an effective 
Hamiltonian, where the kinetic term is replaced by the 
band dispersion — Peierls substitution [lJ-Q . The band 
dispersion is periodic so that the direction of the veloc- 
ity, that is the gradient of the Hamiltonian with respect 
to momentum, oscillates, and the classical trajectories 
generated by the effective Hamiltonian are periodic. 

The high controllability and weak coupling to the en- 
vironment of cold atom systems makes them ideal for 
observation of Bloch oscillation One of the exper- 

imental applications of cold-atom Bloch oscillations has 
been force measurements, in particular of the accelera- 
tion of gravity [l3l - [l8j . that has recently achieved 10~ 7 
accuracy. This application is based on a precise knowl- 
edge of the periodic potential and a large number of os- 
cillation that enable an accurate measurement of the pe- 
riod. The experiments are carried out with very weakly 
interacting Bose atoms to minimize dephasing. 

A cold-atom wave packet undergoing Bloch oscillations 
in an optical lattice with an external force accumulates 
a dynamical phase during its propagation. The dynam- 
ical phase accumulates much faster than the oscillations 
phase, so that a single Bloch oscillation yields a high- 
precision measurement of the external force. However, 
the measurement of a phase requires an interferometric 
setup. For this purpose we propose to use a narrow one 
dimensional potential barrier as a cold-atom wave packet 
splitter to form an interferometer with a Mach-Zehnder 
topology, shown schematically in Fig. [TJ The auxiliary 
mirrors of optical Mach-Zehnder interferometers are not 
required here, since the trajectories in the two arms of 
the cold atom interferometer are curved, and naturally 
rccombine. 

A tunnel barrier for cold atoms can be implemented by 
a narrow blue-detuned laser beam. Coherent tunneling of 



Bosc-Einstcin condensates through a barrier significantly 
narrower than the spatial extent of the condensate has 
been used in experiments to study condensate dynamics 
in a double- well potential [25H27I ]. Here we propose to 
use the same technique in an optical lattice instead of a 
trap. The width of the barrier, equal to several lattice 
spacings, fits the requirements of the interferometer. 

The proposed setup is therefore as follows: a cold-atom 
wave packet is loaded adiabatically into an optical lattice 
with lattice spacing much shorter than the (spatial) wave 
packet width, so that only a single energy band is occu- 
pied. We assume that the band is non-degenerate and 
does not cross other bands. The atoms are also subject 
to a weak external force along a lattice direction. The 
wave packet is accelerated by an external force and im- 
pinges on a narrow potential barrier that is nearly per- 
pendicular to the force direction; it then splits, the wave 
packet fragments propagate in the two arms, and recom- 
bine in a second collision with the potential barrier. The 
interferometer structure is summarized in Fig. [1] High 
sensitivity is derived from the fast oscillations of the in- 
terferometer outputs as a function of the tilt angle of the 
tunnel barrier. 

The interferometer system relies on three small param- 
eters. The principal parameter is s, the ratio of the po- 
tential energy difference across a lattice unit divided by 
the energy band energy scale. The inequality e < 1 guar- 
antees the validity of the semiclassical approximation and 
provides the necessary scale separation between propaga- 
tion and scattering in the operation of the interferome- 
ter. It is also essential for the interferometric sensitivity 
of the phase, since the accumulated phase difference in 
the two arms is proportional to e . The second small 
parameter is 5, the ratio of the lattice spacing and the 
spatial uncertainty of the cold-atom wave packet. 6 -^l 1 
is necessary for the localization of the wave packet in a 
single band, while momentum localization requires that 
e -C 5. The third small parameter is the lattice spacing 
divided by the width w of the tunnel barrier, w has to 
be much shorter than the wave packet width to prevent 
wave packet distortion during the scattering and much 
longer than a to prevent band transitions. The setup 
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FIG. 1: (Color online) A schematic diagram of the cold atom 
interferometer, consisting of a square optical lattice whose 
axes are shown as the x and y axes, an external force F in 
the x direction, and a tunnel barrier (red) tilted by a small 
angle a with respect to the y axis. The initial wave packet 
indicated at the bottom of the figure (blue) at phase space 
point (qi,pi) splits at time U by scattering into a transmitted 
and reflected sub-wave packets, whose classical reference or- 
bits (<7^ T ' ) ,^ T ') and (<7i T \Pi T ' ) ) (respectively) are indicated in 
blue and green (respectively) , that scatter once more from the 
tunnel barrier at times t-r and tR (respectively) and recom- 
bine. The interference of the transmitted-transmitted (TT) 
and refiected-refiected (RR) wave packets is sketched near the 
top of the figure. Some examples of the expected fraction of 
the atoms in the TT-RR output are shown in Figs. [2] and [3] 

conditions are therefore 

e < 5 < — < 1 (1) 

w 

Good visibility of the interference between the two 
wave packets in the outputs of the interferometer depends 
on phase-space overlap. Time-reversal symmetry guar- 
antees that the transmitted and reflected wave packets 
recombine at the beam splitter with full overlap after a 
complete Bloch oscillation if the barrier is perpendicular 
to the force direction, regardless of the band dispersion; 



it follows by continuity that there is significant overlap 
when the tilt angle a of the beam splitter with respect 
to this direction is small enough. The population in the 
outputs exhibits fast oscillations as a function of a with 
a period of order e, shown in Eq. (|2"2"|) and Figs. [2 and 
[3] The range of a with high visibility is of order ei for 
the optimal choice S ~ , so that the number of visible 
oscillation is large, ~ e~2. In the initial, constant vis- 
ibility range of a, the oscillations are purely sinusoidal 
with a period determined by the initial momentum and 
the band dispersion, as shown in Eq. (|24p. 

The wave packet motion in the arms of the interferom- 
eter is planar, undergoing Bloch oscillations in one direc- 
tion and free motion in the perpendicular direction. A 
ID optical lattice is therefore required, with lattice axis 
along the direction of the external force. Nevertheless, 
the interferometer operates equally well in a 2D lattice, 
and such operation offers additional probes, discussed be- 
low. Since the ID-lattice case is obtained as special case 
of the 2D one, we carry out the analysis for a 2D lattice, 
and present the ID special case for the key results. 

The recombination of the wave packet occurs after a 
single Bloch oscillation. However, the interferometer ge- 
ometry implies that the wave packets propagating in the 
interferometer arms cross the barrier once also during 
the Bloch period. The wave packets in general do not re- 
combine during this intermediate crossing, so it is useless 
for the purpose of interferometry. A practical method of 
avoiding these spurious collisions and subsequent degra- 
dation of the interference pattern is to make the barrier 
potential time-dependent. The barrier potential needs 
to be maintained only during the initial and final col- 
lisions, and can be turned off during the intermediate 
crossing of the barrier, since these are well-separated in 
time. We will assume that this plan is carried out in 
our calculations and, for the sake of completeness, point 
out the modifications that arise if the barrier is time- 
independent. 

Bloch oscillations are a wave phenomenon that is not 
particular to quantum mechanics. In addition to cold 
atoms systems they have been experimentally observed 
and studied in guided light systems [ill [2(| ■ These sys- 
tems are typically paraxial, and since the paraxial wave 
equation is identical to the Schrodinger equation with the 
propagation direction playing the role of the time coor- 
dinate, the propagation of a quantum wave packet in a 
periodic potential and of a beam in wave guide array are 
described by the same analysis. This result also holds 
for the present study; for concreteness we express all our 
results in cold atom terms, as it is the more natural ap- 
plication. 

In addition to the spatial interferometry considered 
here, waves propagating in a periodic medium can split 
and interfere in the band degree of freedom [21M241 ] . In- 
stead of a tunnel barrier, this requires a coupling between 
bands, possibly by Landau-Zener tunneling. This type of 
interference does not directly measure the external force 
unless the periodic structure is controlled. 
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The interferometer operation is based on two building 
blocks: Semiclassical wave packet propagation [29| in the 
arms of the interferometer, that is pre sented in Sec. [HI 
and quantum normal form analysis [301 ] of the scattering 
of the wave packets off the tunnel barrier, that is the sub- 
ject of Sec. Mil The results of the wave packet propaga- 
tion analysis and wave packet scattering are combined in 
Scc. lIVI to yield expressions for the interferometer output 
as a function of the tilt angle. In Sec. [V] we discuss ap- 
plications of the interferometer force measurement and 
as a probe of the optical lattice structure. It is shown 
that the interferometer is experimentally realizable for a 
standard setup with gravity as the external force 



II. WAVE PACKET BLOCH OSCILLATIONS 

The propagation of the wave packets in the arms of 
the interferometer is governed by the combination of the 
optical lattice potential, and the external force F directed 
along a lattice axis, that we choose as the x axis, so that 
the single-atom Hamiltonian is 

^ + W£> (2) 

where q and p are the atom's position and momentum, 
respectively, and V per is the lattice potential. We assume 
that the inter-atom interactions are negligible, and that 
the many-atom state is a product of identical one-atom 
states 

The optical lattice is either one- or two-dimensional. In 
either case, we suppose that the wave packet is fully sup- 
ported in one of the energy bands of the (unforced) lat- 
tice potential whose energy dispersion is £(p x ,p y ), where 
p is the quasi-momcntum, and measure energy in units 
of U, the energy band width. If the optical lattice is one- 
dimensional then £(p x ,p y ) = £i(p x ) + 2^P«) where £\ 
is the one-dimensional dispersion, and p y is the y com- 
ponent of the ordinary momentum. We assume that the 
energy band is nondegenerate and separated from the 
other bands by a gap of order U. We also choose -, 
where a is the lattice spacing, as the unit of momentum, 
so that £ is periodic with period one in p x . The exter- 
nal force is weak, in the sense that e = ^y- -C 1. In 
the following we will write simply momentum instead of 
quasi-momentum, when there is no risk of confusion. 

The initial single- atom state \i) is a wave packet with 
mean position qi and momentum pi and momentum un- 
certainty S <C 1. It follows that the position uncer- 
tainty is much larger than a. In explicit calculations 
we let \i) be a Gaussian wave packet with wavefunction 

(pit) = cexp(-^=f^- -iti -p). 

Under these conditions, an effective classical dynamics 
can be applied as a systematic approximation for the 
propagation of the wave packet. An effective Hamil- 
tonian H e e determines the in-band dynamics, derived 
from a semiclassical expansion in e, whose leading term, 



Peierls substitution pj, is obtained by replacing the ki- 
netic term in the Hamiltonian with the band dispersion, 
so that 

H cS = £{p) - q- F (3) 

and letting q and p be canonically conjugate [q n ,Pm] = 
ihS nm . q is a coarse-grained position observables with 
an uncertainty much larger than a. We assume that the 
optical lattice is space-reflection invariant in x and y, so 
that the sub-leading term in H e s is 0(e 2 ) and therefore 
negligible. 

The semiclassical propagation of the wave packet gen- 
erated by H e s consists of three elements, a phase-space 
drift along the classical trajectory of the center of the 
wave packet (qt,Pt), a deformation (squeezing) of the 
wave packet by the linearized flow at the wave packet 
center, and an overall phase accumulation proportional 
to the classical action of the wave packet center trajec- 
tory [28l.l29l]. Each part is realized as a unitary operator, 
so that the wave packet at time t f can be represented by 

\f) = e^T(q f ,p f )M(S fl )THq ti m) (4) 

where 

1 f f 

7=--y (±(p-dq-q-dp)-H cS dt) 

is the dynamical phase, 

T(q ,Po) = exp(±|(g-p -f- 3o)) 

is the Hciscnbcrg phase-space shift operator, and M{S) 
is the squeezing (metaplectic) operator derived from the 
symplcctic phase-space deformation matrix 




with block components S qq = dqf in /dqi t7n , S pq = 
dpf t n/dq it m, S qp = dqf tn /dpi, m , and S pp = dpf tn /dp itTn . 

The problem is therefore reduced to the calculation of 
the classical trajectories of H c ^. The effective Hamilto- 
nian is integrable, and its flow can be calculated explic- 
itly. The Hamilton equations are 

8 t q = V£(p) (5) 
d t p = F (6) 

The momentum equation implies that pt — Pi + Ft: 
the y position equation q t . y = qi. y + J t . dt'v y (t'), where 
v y (t) = d y £(pt), and conservation of energy implies that 
qt,x = q%.x + The periodicity of £ im- 

plies that the wave packet performs Bloch oscillations in 
the x direction with period y , while moving freely in the 
y direction. 

It follows that the S qq and S pp blocks of the deforma- 
tion matrix are unit matrices, the S pq block is the zero 
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matrix, and the elements of the S qp block are 
dgf, x 1 



dpi 



F 



(vx(pf) - v x (pi)) 



(7) 



d Pl , y F 



— / dp x dyV y {p Xl P l .y) 



(9) 



The expression for the dynamical phase can be simplified 
to 



1 



U 



f 



W = 2^(P> • Qf ~Pi ■ Q%) ~ ftp I dp x £(p x ,p hy ) (10) 



In our choice of units the coefficient -re — — , so that 

nF e 1 

the classical limit, where h — > keeping all classical quan- 
tities fixed, is equivalent here to the limit e — > 0, and e is 
the effective Planck constant. In this limit the dynami- 
cal phase diverges, giving rise to the high interferometric 
sensitivity, as the small wavelength does in classical op- 
tical interferometry. 



III. WAVE PACKET SPLITTING OFF A 
TUNNEL BARRIER 

The atom wave packet is split by tunneling through a 
one-dimensional potential barrier or a potential trench, 
that is, a modulation of the optical potential of the form 
Ub(q ■ h), where h is the unit normal to the beam splitter, 
and Ub is localized on a scale w of magnitude a <C w <C 
a/5. 

For the purpose of analyzing the scattering of the 
wave packet by the beam splitter, we ignore temporar- 
ily the external force, and consider a standard scatter- 
ing problem of an incoming unit amplitude plane wave 
with quasi-momentum p in an energy band with disper- 
sion £ , that splits into a transmitted wave with quasi- 
momentum p with amplitude t and a reflected wave with 
quasi-momentum ff R ^ with amplitude r. The inequal- 
ity a <C w guarantees that the scattered wave remains 
in the same energy band. Lattice translation symmetry 
implies that tangential component of quasi-momentum 
is conserved, pt = p\ , and energy conservation implies 
that £(p) = £(jf R ^); we assume that the band dispersion 
is such that there is a single solution p^ RS> , other than p, 
to these two conditions. 

The initial quasi-momentum and wave packet splitter 
orientation define a classical band of allowed energies be- 
tween £_ = min p?i £(p+p n h) and £+ = max Pn £(p+p n n). 
The condition a <C w also implies that unless £(p) — £_ 
is close to max Mb or £{p) — £+ is close to minwf, ei- 
ther \t\ <C 1 or | j- 1 <C 1. For concreteness we choose 
to study scattering from a potential barrier and let 
\£(p) — £- — max it;, | C 1 so that both t and r are appre- 
ciable to obtain good interferometric visibility. We also 



assume that the initial momentum is such that £(pi) is 
far from the edges £ ± of the interval of allowed energies. 

The typical barrier maximum and dispersion minimum 
are quadratic, and under the conditions laid out above 
the scattering amplitudes are determined by local behav- 
ior of £ and Ub near their extrema [30l | . Denoting by p m 
and q m the abscissas of the extrema, we define the clas- 
sical action 



I = 



£{p)- £{p m ) - Ubiq^ 



-dl£(p m )9Mq 



(ii) 



and the standard theory of scattering from a parabolic 
barrier gives [U 



t «=^^- r (2-n) e " 

e ii(logi-l) x J 

i 1= r - — i— e 



(12) 
(13) 



A scattering wave packet is a superposition of sta- 
tionary scattering states \i) = J dpip{p)\p). After the 
center of the incoming wave packet collides with the 
barrier it splits into a transmitted wave packet \i) T = 
/ dpt(p)i/j(p)\p) and a reflected wave packet \i) R = 

f dpr(p)^(p)y/\fetT\\p(V), where J = || <R) is the Ja- 
cobian matrix of the reflection transformation. Since the 
incoming wave packet is concentrated near pi, the trans- 
mitted wave packet is approximately 



dpt(pi)i>{p)\p) = t(pi)\i) 



(14) 



but the dependence of the scattering amplitudes on p 
causes a position-space shift and distortion of the wave 
packet. Nevertheless, the shift and distortion are weak if 
the t and r do not change appreciably over the range of 
momenta that support the wave packet, or equivalcntly 
if the barrier width w is much smaller than the posi- 
tion uncertainty j (32|, as we assume. In this case the 
position and momentum shifts incurred by the scatter- 
ing are of 0(w) and 0(^5 2 ) (respectively), negligible in 
comparison with the respective uncertainties; the relative 
deformation is even smaller, of 0((^) 2 ). 

By the same reasoning the reflected wave packet can be 
written as \i) R = r(pi) J dpt[)(p)^\ det J\\p^ R '). Since the 
momentum integration is localized, J can be evaluated at 
Pi, and the reflected momentum is approximately piR) = 
jf^ + Ji{p — Pi), where Ji = J{pi). The reflected wave 
packet is therefore 



\i) R = r(p i )T(fi)T(p<V)M(Mr 1 o ) Tt(pl)T t ((f . )H) 

(15) 

Finally we reconsider the effect of external force. In 
this case the quasi-momentum is no longer a good quan- 
tum number, so we label the stationary states by the 
value of p at the barrier. These states are propagating in 
a large interval of O(-) size around the barrier and can 
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therefore be used as a basis for a scattering theory for the 
wave packets considered here, with position uncertainty 
of 0(f) «C | . The preceding arguments and Eqs. ([T2T - 
I15[) arc valid also with the weak uniform external force 
up to small corrections. 



IV. ANGLE INTERFEROMETRY 

The wave packet interferometry takes place in three 
steps. An initial wave packet |z) of mean position qi 
and momentum pi is first split by scattering from the 
wave packet splitter into a transmitted wave packet \i) T 
and a reflected wave packet with mean position q\ 

-L R) 

and momentum p\ . The transmitted wave packet then 
undergoes Bloch oscillation in the x direction and prop- 
agates freely in the y direction before impinging on the 
beam splitter again at time tx and position with mean 

momentum pi . The wave packet \c) T then splits again 
into a transmitted-transmitted (TT) wavepacket \c) TT at 
qi T and p^ T ^ and a transmitted-reflected (TR) wave 

packet \c) TR at q^ R ^ andpi T . The states and variables 
related to the propagation and splitting of the reflected 
wave packet arc similarly defined, sec Fig. [TJ 

The system has a Mach-Zehnder geometry with two 
outputs, one interfering TT with RR wave packets, and 
the other interfering TR with RT wave packets. The in- 
terferometer is sensitive to the difference between the dy- 
namical phases accumulated by the transmitted and re- 
flected wavepackets during their propagation. This phase 
difference changes fast as a function of the tilt angle a 
between the beam splitter and the y axis. 



A. Wave packet interferometry 

We now calculate the interference pattern, concentrat- 
ing from this point on the TT-RR output. The inter- 
ference takes place after both the transmitted and the 
reflected wave packets have collided again with the beam 
splitter. When the tilt angle a is positive > tx, so 
that tf = 4r. The final state of the TT wave packet is 
therefore 

x A/(4fV(<?1,K)N} (16) 
with momentum-space wave function 



JTT) 



i r-lTT) -/TT) _ _ , i -ITT) - 

X pn^f 'Pi -3*'P0 e -*9/ P 



1 f-* -i TT )\ n(TT\ t- ATT)-, 



for C*( TT ' = 1 + i^-S^JJ, and the final state of the RR 
wave packet is 

|/>*« = | c )*a = ri^yip^n^) 

xT(e\^)M(Sif)T^,^) 

x T(<M)T(ti R) )M{ Mr 1 I )Ti(jH)T\q})\i) (18) 

with wave function 



-(H) 



dct(J I )dct(4 i?) ) 

(p-p< rr ) ) . c' hr ' ■ (p-pi RR} ) 



x e « 



(19) 



for 



= ((^)*)- 1 ((j/)- 1 jt 1 + i 2 4s^ q )(A R) )- 1 - 

Conservation of energy and transverse momentum imply 
that in Cartesian coordinates 



J, = ( V v *,v 
sin a cos a 



j(R) = VI 



(RR) (RR) 



— sin a cos a 



"\>* l ' v I (20) 

- sm a cos a / 

„(H) ,.(«) , 

Uc.x u c,y \ (21) 

— sin a cos a } 



The state in the TT-RR output is If) 11 + \c) HH , so 
that the probability that an atom is detected there is 
\t(/PWi)\ 2 + \r(^ c R) )r(^)\ 2 + 2Re TT (f\c) RR . The in- 
tcrfcrcnce term takes the form 

rr (/|c> M = t(^rt(ft)'r^)r(ft) 
x e 

- ! rj <TT) ,?< TT )_o. e-1 J.R) - -OR), 

e l 1 '/ p / 9 » p »> e f -Pc ~H-p\ ) 

v /dot(J 1 )dct(J^ R) )dct(i((C< TT ))*+C< RR ))) 

1 (/ 

x e 



r) .(C (T-r ')*-p?" r) +p< HR) .c< RR ) -p< RJi> ) 



(22) 



where 



(17) 



The phase of the interference term changes as a function 
of a on a scale of e, while the amplitude and the frequency 
of phase oscillations change more slowly. These proper- 
ties are visible in the interference patterns shown by thick 
blue lines in Figs. [5] and [3] for one- and two-dimensional 
optical lattices, respectively. The interference patterns 
were calculated for two values of e, 10~ 3 and 10 -4 . 
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FIG. 2: (Color online) Probability of observing an atom in the TT-RR output of the interferometer as a function of the tilt 
angle a. The atoms propagate in a one-dimensional tight-binding potential period-two super-lattice with hopping energy ^U, 

on-site potential ±(7, and lattice spacing a — 2nh)J~^. The thick blue line was calculated with the full expression (I22|l . and 
the thin violet is the small angle approximation (|24|) . The initial momentum space width 5 — \\fe and the barrier width 
w = 2tt"^?" The left and right panels show the interference patterns for e = 10~ 3 and e = 10 -4 , respectively. 
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FIG. 3: (Color online) The interferometer output for propagation in a separable two-dimensional optical lattice that is a 
superposition of x and y potential super-lattices having the same properties as the one used in Fig. [2] The significance of the 
axes, the curves, the initial state, the width of the barrier, and the values of e are the same as in Fig. [2] 



The results presented in this subsection were derived 
under the assumption that the T and R wave packets do 
not scatter before the recombination event. Since the tra- 
jectories of the wave packets in the interferometer arms 
cross the barrier line this requires temporal control of the 
barrier amplitude so that it vanishes at the intermediate 
crossing times of the T and R wave packets, and kept 
constant at the initial splitting time and final recombi- 
nation times. If the barrier amplitude is not controlled 
temporally, further splitting would make the final state 
a superposition of six (rather than four) wave packets 
or more. Nevertheless the interference in the TT-RR 
and TR-RT channels persist with the modification that 
the state |/) TT is multiplied by an additional transmis- 
sion amplitude t(p\ ) compared with Eqs. (fTo) and (JT7J), 

-IT) 

where p\ is the momentum of the T wave packet when 
it crosses the barrier, and similarly \c) RR gains an ad- 
ditional factor of t(px) with respect to Eqs. (fT8)l and 
(TK)|) . In the following we assume that the barrier poten- 



tial is controlled appropriately and therefore omit these 
additional factors that reduce the amplitude of the inter- 
fcrometric oscillations in Eq. ([2^|) . but otherwise do not 
affect the results. 



B. Interferometry for small angles 

The interference is effective when the interference term 
is large, and this requires that the wave packets 
|/) TT and \c) RR overlap in position as well as in mo- 
mentum. For most tilt angles q^" 1 '^ ^ qi RR ^ and p"J T ^ ^ 
pc ■ However, if the beam splitter is aligned with the 



(ft) 



-Pi,x, P 



(ft) 



Pi 



and 



y axis then by symmetry p, 
v y(—Px,Py) = Vy{p x ,p y ) 1 so that after one Bloch period 
the two sub-wave packets meet at the beam splitter at 

(T) (ft) (T) (ft) ilnc/ \ v, 

q c , x = q c J = Qi, x , q ft y = q f , y = 1i,y + FOy£{Pi,y) where 
8{pi. y ) = J Q dp x v y (p x ,pi^y). There is therefore full over- 
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lap between the TT and the RR wave packets for a = 0, 
and high interferometric visibility for small enough a. 

Since a is small, it is possible to approximate the classi- 
cal data in the interference term by its Taylor expansion. 
The 0(a) approximation gives 

TT (f\c) RR = (t(jS<)V(ft)) a (l + a(ht + h 2 )) (23) 
where 

h 2 = -i9 Q (dct(J 4 )dct(ji i?) )dct(i((C( TT ))* + C( fl ^)) 

The second correction to the interference term ah 2 
is 0(a) and therefore negligible for small a; similarly, 
a/ii = O(^-a) is also negligible throughout the region of 
validity of Eq. (f2"3"|) , where, therefore, 

(/|c) flil = {fr) 2 e s ^ BB ^ a (24) 

using the shorthand t = t(pi), r = r{pi). It follows that 
for small enough angles the probability of observing the 
atom in the TT-RR output oscillates sinusoidally as a 
function of a on a fast scale of 0(e) and large visibility. 
The thin violet lines in Figs. [5] and [3] show the interference 
pattern with the approximation Eq. (|24p that becomes in- 
creasingly accurate for smaller a. If the band dispersion 
is known, these oscillations can be used to measure the 
force F with high precision, as discussed further below. 

The domain of validity of Eqs. and (glj) is limited 
by the size of the neglected 0(a 2 ) terms. These terms 

2 2 2 c2 

arise in the combinations and 9L p- in Eq. (I22j) . 

In particular, the first and third combinations determine 
the scale of decay of visibility as a result of momentum 
mismatch and position mismatch between the TT and 
RR wave packets, respectively. It follows that the domain 
of validity is 

a < min(5,£5, f) (25) 

Evidently, the largest range is obtained when 8 ~ , 
that is when the momentum and position uncertainties 
are balanced, so that the domain of validity is a <C e 3 . 
In this case, since the oscillation period is O(e), the num- 
ber of equal period oscillations up to a fixed tolerance is 
0(e-2). When a is comparable with e 2 the visibility of 
the interference pattern decreases and the oscillation pe- 
riod changes, as shown in Fig. [5] and [21 that also demon- 
strate that the number of high- visibility constant-period 
oscillations increases when e becomes smaller. 



V. CONCLUSIONS 

An experiment of angle interferometry would detect 
the fraction of atoms scattered in the TT-RR arm for 
variable tilt angle of the tunnel barrier, keeping the rest 



of the parameters fixed. Writing Eq. ([2"4")) in physical 
units, this fraction would be approximately equal to the 
probability 

(f\c) RR = |i| 4 + |r| 4 + 2Re((rr) 2 e 2j 7tF^(p^)) (26) 

With the optimal choice of initial uncertainty this ex- 
pression is valid for values of a significantly smaller than 
e~ ' . A time of flight measurement can yield the squared 
magnitude of the full momentum-space wave function, 
the TT-RR part of which is given by Eqs. (fTT)l and (fT§jl . 

The interferometer period can be further simplified 
when the optical lattice is separable, a case that includes 
one-dimensional lattices, obtaining 

(f\c) RR = |i| 4 + M 4 + 2 Re((t*r) 2 e 2l ^ a ) (27) 

where Vi^ y is the constant velocity in the direction per- 
pendicular to the force. 

If gravity is used as the external force acting on atoms 
with mass m in lattice with spacing a then 



where g is the acceleration of gravity, and u is the band 
width in units of the recoil energy . It follows that 
the interferometric condition e <C 1 is experimentally ac- 
cessible. For example, for Sodium atoms in a ^589 nm 
lattice e = -3.3 x 10 -2 ; a lattice depth of 0.25 recoil en- 
ergy, that can be set up with sub-mW lasers (HI, gives 
u w 0.4, and e = 8.4 x 10~ 2 . Spatial localization of 3 /im 
and barrier width of 1 /xm can be chosen to satisfy the 
basic inequalities Eq. ([1]). 

A natural application of the interferometer is force 
measurement, that can reach high precision thanks to 
the high sensitivity of the interference pattern, and for 
this application a one-dimensional lattice suffices, and 
the force can be deduced from the interferometer output 
using Eq. (|26|) . The statistical error in a single measure- 
ment is of the order of the inverse square root of the 
number of atoms, that can be 10 -3 in typical experi- 
ments. If the measurement is repeated for one hundred 
values of a the statistical error is reduced by another fac- 
tor of 10. The induced measurement error in the force 
is reduced by the number of visible interferometric os- 
cillations that is of order y/e, giving an error estimate 
of 10~ 5 . The interferometer measures the combination 
rather than the force itself, so that high-precision 
measurement of the transverse speed is needed to reach 
this value of force measurement sensitivity. This sensi- 
tivity is not as good as that of existing interferometers 
that rely on the classical Bloch oscillation phase, but the 
dynamical phase interferometer has the advantage of re- 
quiring a single Bloch oscillation, facilitating much faster 
measurements. 

The interferometer can also be used to characterize the 
dispersion of two-dimensional lattices with a given exter- 
nal force. This application can become especially inter- 
esting if the lattice breaks space-reflection symmetry so 
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that the interferometer can measure a non-trivial Berry 
phase [HI. This question, however, requires a more ac- 
curate analysis that is beyond the scope of the present 
paper. 
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